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Abstract 

Our goal is to find closed form analytic expressions for the solitary waves of nonlinear 
nonintegrable partial differential equations. The suitable methods, which can only be nonper- 
turbative, are classified in two classes. 

In the first class, which includes the well known so-called truncation methods, one a priori 
assumes a given class of expressions (polynomials, etc) for the unknown solution; the involved 
work can easily be done by hand but all solutions outside the given class are surely missed. 

In the second class, instead of searching an expression for the solution, one builds an 
intermediate, equivalent information, namely the first order autonomous ODE satisfied by the 
solitary wave; in principle, no solution can be missed, but the involved work requires computer 
algebra. 

We present the application to the cubic and quintic complex one- dimensional Ginzburg- 
Landau equations, and to the Kuramoto-Sivashinsky equation. 

Keywords: solitary waves, complex one-dimensional Ginzburg-Landau equation, Kuramoto-Sivashinsky 
equation, complex Swift-Hohenberg equation, Briot and Bouquet equations, elliptic function, 
genus, Painleve property, meromorphy, truncation. 
PACS 1995 : 02.30.-f 05.45. +b 42.65.-k 47.27.-i 
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1 Introduction 



Many nonlinear partial differential equations (PDEs) encountered in physics are autonomous, i.e. 
do not depend explicitly on the independent variables x (space) and t (time). In such a case, 
they admit a reduction, called traveling wave reduction, to an autonomous nonlinear ordinary 
differential equation (ODE), defined in the simplest case by u(x,t) = £/(£),£ = x — ct, with c a 
constant speed. A solitary wave is then defined as any solution to this nonlinear autonomous ODE. 
Physically relevant solitary waves must satisfy some decaying condition when £ goes to ±00. 

The distinctive feature of this chapter is to explain the methods to find closed form expressions 
to these solitary waves when the PDE and the reduced ODE are algebraic and nonintegrable. 
These solitary waves may have the topology of a front (for instance tanh), a pulse (for instance 
sech) , a source, a sink, etc, but we will not discard an apparently physically uninteresting solution, 
because it might appear interesting to another field. 

Why "algebraic"? This only excludes equations impossible to convert to an algebraic form. For 
instance, the sine-Gordon equation u xt — sinu = is not excluded because it is algebraic in e m . 

Why "nonintegrable"? Because the integrable ones (nonlinear Schrodinger (NLS), coupled NLS 
in the Manakov case, etc) are "easy" to solve using powerful tools like the inverse spectral transform 
(1ST) 1 . The difficulty with the nonintegrable equations is the absence of a general method to 
achieve the goal. 

Why "autonomous" ? Irrelevant for the truncation methods (section EJ , this restriction is 
essential for the mathematical method of section[7J Physically, this is not an important restriction, 
since many interesting PDEs are autonomous, see the examples below. 

Why "closed form expressions"? Because a solution represented by a series can be misleading 
(illusoire, used to say Painleve). Consider for instance a chaotic deterministic dynamical system 
for which no analytic solution exists. Around a regular point, it admits a solution represented by a 
Taylor series, but one can conclude nothing before some analytic continuation has been performed. 
On the contrary, the Laurent series around a movable singularity (i.e. one whose location depends 
on the initial conditions) provides some constructive information (see section lXTJl about the (global) 
integrability of the equation. 

The methods described here are all based on the a priori singularities [23 of the solutions of 
the given ODE. In particular, we do not consider the group theoretical methods |4()j . 

These methods can be applied mainly to dissipative equations of importance in physics, non- 
linear optics, mechanics, etc. Our specific examples are the following. 

1. The one-dimensional cubic complex Ginzburg-Landau equation (CGL3) 

iA t + P A xx + q\A\ 2 A~i 1 A = 0, pqj^O, Tm(p/q) ^ 0, (A,p,q)eC, 7 eK, (1) 

(and its complex conjugate, i.e. a total differential order four), in which p, q, 7 are constants, 
a generic equation which describes many physical phenomena, such as the propagation of a 
signal in an optical fiber [2J, spatiotemporal intermittency in spatially extended dissipative 
systems We will restrict ourselves to the CGL3 case properly said lm(p/q) 7^ 0. 

For analytic results on two coupled CGL3 equations, see 

2. The Kuramoto and Sivashinsky (KS) equation, 

ft + vfxxxx + of xxx + Vfxx + ffx = 0, f eC, (v, b, n) G ft, V ^ 0. (2) 

in which f,b,fi are constants. This PDE is obeyed by the variable tp = aigA of the above 
field A of CGL3 under some limit |321 12(J| . hence its name of phase turbulence equation. 

3. The one-dimensional quintic complex Ginzburg-Landau equation (CGL5), 

iA t + P A xx + q\A\ 2 A + r\A\ 4 A-ijA = 0, pr^O, Im(p/r) ^ 0, (A,p,q,r)eC, 7 eft. (3) 

4. The Swift-Hohenberg equation [371121] 

1 A t + bA xxxx + P A XX + q\A\ 2 A + r\A\ 4 A - ijA = 0, br ^ 0, (A, b,p, q,r,)eC, 7 G ft, (4) 
in which b,p, q, r, 7 are constants. 
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For the CGL3, KS, and Swift- Hohcnbcrg equations (with one exception, KS with b 2 — 16/i^), 
all the solitary wave solutions \ A\ 2 = /(£), ip = $(£ ), £ = x — ct, which are known hitherto are just 
polynomials in tanh 

So, two natural questions arise: 

1. Is it possible or impossible that other solitary waves exist? 

2. If such other solutions may exist, can one find not just a few more but all of them? 
Let us from now on denote the reduced ODE as 

..,u»=0, ' = ±,Z = x -ct. (5) 

and let us assume that it is also nonintegrable. 
The chapter is organized as follows. 

In section [2 one recalls the analytic expressions of the known solitary waves of the examples. 
This list, to be retrieved or augmented by the singularity based methods, will allow us to rate the 
efficiency of the various methods. 

In section one investigates the amount of integrability of the equation, by applying the so- 
called Painleve test. More specifically, one checks the existence of particular solutions which admit 
a local representation as a Laurent series. This allows us to count the gap, strictly positive because 
of the assumed nonintegrability, between the differential order of the ODE and the maximal number 
of available integration constants. 

In section^] we discuss the choice of the suitable dependent variable to be used in the subsequent 
sections. 

In section|Sl one tries to obtain a global representation for the local information (Laurent series) 
previously found. One introduces the distinction between two main classes of methods, according 
to the following criteria: i) computations easy enough to be carried out by hand, ii) generality or 
particularity of the expected solution. 

The sectionals devoted to the first class of methods, which are known as "truncation methods" . 
The input is a class of a priori expressions for u (usually polynomials), in some intermediate variable 
X which satisfies a given first order ODE (e.g. Riccati, Weierstrass, Jacobi). Then, by a direct 
computation, easy to carry out by hand, one checks whether there indeed exist solutions in the 
given class. The solutions with a simple profile (such as tanh for a front, sech for a pulse), are 
easily found by this class of methods. 

In the second class of methods presented in section rather than directly looking for the 
unknown solution 

u = /(£-&), (6) 

in which £o is an arbitrary complex constant, one looks as an intermediate information for the first 
order nonlinear ODE 

F(u,u')=0, (7) 

obtained by eliminating £o between Q and its derivative, in which F is as unknown as /. Indeed, 
provided that / is singlevalued, by a classical theorem recalled in Appendix, there is equivalence 
betwen the knowledge of the solution / and that of the subequation F which it satisfies. The way 
to obtain the subequation F is to require that it be satisfied by the Laurent scries obtained in a 
previous step. 

The difference between the two classes of methods is the following. The solutions found by 
the first class of methods can only be a subset of those found by the second class. However, the 
computations involved can easily be performed by hand for the first class, while for the second 
class a computer algebra package is highly recommended. 

2 The known solutions of the examples 

None of the expressions listed below represents the largest analytic solution which one could find, 
and their distance to this largest, yet unknown, analytic solution will be computed precisely in 
section |21 
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2.1 CGL3 

The traveling wave reduction of Q 

A(x,t) = y/M^)e i( -~ wt + £ = x-ct, (c, u, M, ip) G TZ, 



M' 1 



2M AM 2 



M j.- » > 2 m ■ < M ' ■ c M ' j. 1 ( >m u q at ^ n 

+ Z<^ -tp + up — - t - TJ + 'iCif +LU) + -M =0, 

M 2p M p p p 



(8) 
(9) 



introduces two additional real constants (c, uj) and it is convenient to define the six real parameters 



7 I • i 9 • 1 . . 7 + » . 1 22 
d r + id; = ~, s r - iSi = -, ,g. r + = — h -c s r . 



P 



P 



P 



(10) 



In the CGL3 case properly said di 7^ to which we restrict here, only three solutions are 
currently known. Denoting Aq and a two real constants defined by the complex equation 



(-1 + ia)(-2 + ia)p + Afo 
these three solutions are the following. 

1. A heteroclinic source or propagating hole 0] 



0. 



A = A Q 

ij — LU 
P 



iqpi 



k , ife 

— tanh — £ 

2 2" 2(l-ta)p\p\ 2 d i 
k 2 



i[a Log cosh — £ 



5 — c£ — uit\ 

2\p\ 2 d % 



(12) 



= 1^) -(2-3ia)-, 



in which the velocity c is arbitrary. Indeed, the real and imaginary parts of the last equation 
define the value of u> and a linear relation between c 2 and k 2 , see jSJ Eq. (79)]. 



2. A homoclinic pulse or solitary wave [3D] 

A = A (-ifc sech kx)e l \- a Lo S cosh fca; ~ ^ , 
!l^ = (l- ia ) 2 fc 2 , c = 0. 



(13) 



3. A heteroclinic front or shock 



An 



17 — UJ 

p 



tanh — £ + e 



i[a Log cosh — £ 



3gr + qpj 
6b| 2 



e z = h 



(14) 



A; 2 fc PiC 

T» 2 =£ 



6|p| 



2 - 



None of these three solutions requires any constraint on p, 7, and they depend on an additional 
sign resulting from the resolution of Ijll)). 



2 _ 3(3d r +£iA) 



n - 3dr + £lA , A = ,/9d 2 + 8rf 2 , e 2 = 1. 



(15) 



2d 2 2d, 

In all of them M is a much simpler expression, namely a second degree polynomial in 

r = (fc/2)tanhfc£/2, k 2 G ft. (16) 

Therefore, if one wants to extend the three above solutions, it is advisable to eliminate tp between 
the system of two real equations equivalent to J^J, 



( M" M' 2 



2M 4M 2 

f +tp w- Sr 



tp - Sj 

( cM 1 



+ 7 + s r (ctp' + u) + d r M = 0, 



\2M 



f cM' 
\2M 

+ 7] - Si (op' + w) + diM = 0, 



(17) 
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which results in 



cs r G'-2csiG / , cs r \ 2 G 

* =— + 2M2(g r -d t My ^~—) = W (18) 

(C - 2c Sl G) 2 - AGM 2 {d l M - g r f = 0, (19) 

1 1 cs- 

G = -MM" - -M' 2 - -^MM' + d r M 3 + g t M 2 , (20) 



and to concentrate on the single third order equation (|19fl for M = \A\* 
2.2 KS 

The traveling wave reduction is defined as 

,2-1' 



ip(x, t) = c + «(£), £ = x — ct, 
which integrates once as 



bu" + ptv! + — 



0, (y,b, rieTl, u^O, (21) 



u 2 

vv!" + bu" + fxu' + — + A = 0, (22) 

in which A in an integration constant. It has a chaotic behavior |22j . and it depends on two 
dimcnsionless parameters, b 2 /{fiv) and vAj /x 3 . 

The known solutions are one elliptic solution, six trigonometric solutions, and one rational 
solution. 

The unique known elliptic solution exists for one constraint between the parameters is, b, [i of 

the PDE d3 ng, 



bu a 2 13/i 3 + vA 

b 2 - 16^ = 0: u = -60vp'-15bp-£, 92 = ^, 9s = ^ ■ 

in which p is the elliptic function of Weierstrass, defined by the ODE 

p /2 = Ap 3 - g 2 p - g 3 . (24) 

The six trigonometric solutions all of them rational in e fc ^, exist at the price of one 

constraint between is, b, /i and another one on A, 

u = ^-IS^+f^-SO^-^V + ^-^^ + T^' ^ 

\19 4x19^/ 2 32 x 19i/ 2 4x19^ 

fc k 
t = -tanh -(£-&), 

the allowed values being listed in Tabled 



Table 1: The six known trigonometric solutions of KS, Eq. I|22|) . They all have the form l|25() . 
The last line is a degeneracy of the elliptic solution (|23|) . 



&7(H 


vAj '/i 3 


vk 2 / n 





-4950/19 3 , 450/19 3 


11/19, -1/19 


144/47 


-1800/47 3 


1/47 


256/73 


-4050/73 3 


1/73 


16 


-18, -8 


1, -1 



Finally, the unique known rational solution 

b = 0, (a = 0, A = u = I20v{£ - ^y 3 , 

6 



(26) 



is a limit of all the above solutions. 

A nice property common to all those solutions is to admit the representation 



V Log ip + constant, V = 60v 



156 



d 2 15(16^ -b 2 ) d 



d£ 3 d£ 



7Qv df' 



(27) 



in which i\) is an entire function (i.e. one without any singularity at a finite distance) whose ODE 
is easy to build, respectively, 

(- Log V)'" 2 - 4(- Log ^)" 3 + .92 (- Log if,)" + .93 = 0, 



k 2 



(28) 
(29) 

ip" = 0. (30) 
This linear operator 2?, which captures the singularity structure, is called the singular part operator. 

2.3 CGL5 

The traveling wave reduction is quite similar to that of CGL3, so we do not repeat it. Again, Aq 
and a denote two real constants defined by the complex equation 



(-1/2 + ta)(-3/2 + ia)p + A%r = 0, 



and the convenient constants arc 



r 1 7 + iw 1,, 

e r + id — -, s r - tSi = -, g r + iQi = h -c s . 

p p p 4 



(31) 



(32) 



In the CGL5 case properly said e, ^ to which we restrict here, only two solutions are currently 
known. 

1. A heteroclinic front or shock |35) . 

f A = A ((fc/2)(tanhfc£/2 + e)) 1/2 e^ a Lo S cosh ^/ 2 + ^ ~ £ 2 = 1, 



(-1/2 + ia) [i(c - 2pK ) + 2e(-2 + ia)pk]p~ A 2 a q = 0, 

2 / x 2 



17 — cj 
P 



K h e(l - ia)Jfe/2 

2p 



(33) 



The number of constraints among (p, 9,r, 7) is either two (case Sj = 0, c arbitrary), or zero 
(case Si 7^ 0, with a fixed velocity). 



2. A homoclinic source or sink 

, , k sinh ka 

A = A 



''n 



1/2 , 



[a Log(cosh fc£ + cosh ka) + K£ — uit] 



cosh fc£ + cosh ka 

c — 2pK — 0, which implies cpi = 0, 

(-1/2 + ia) [-2fcp (l - ia) + 2(-2 + ia)r ] p - ^9 = 0, 



(34) 



27 — w 
P 



+ (1/2 - iafk 2 + (3 - lOia - Aa 2 )kn r + (3 - 8ia - 2a 2 )r 2 /2, 



I r (r 2 + 2fc^ r + fc 2 ) = 0. 

in which K, fc 2 , ro, Mo — cothfca are real constants. The number of constraints among 
(fi 9; r ; 7) is either two (case Sj = 0, with c arbitrary), or one (case s, 7^ 0, c = 0). 



Each of these solutions depends on two additional signs arising from the resolution of (|31H 

Al=e 2 



»2e r + eiA 2e r + £iA r~T ~~ 7 , 2 

-, a = , A = x /4e2 + 3ef, e\ = e\ = 1. 



2e, : 



(35) 
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2.4 Swift-Hohenberg 

Again, Aq and a denote two real constants denned by 

(-1 + ia)(-2 + ia)(-3 + ia)(-A + ia)b + A^r = 0. 



(36) 



In the case Im(r/&) ^ to which we restrict ourselves, only two solutions seem to be currently 
known. 

1. A stationary front gl Eq. (127)] 



' A = A (fe/2)tanhfcr/2 e *[« Log cosh fcx/2 - wt] ] 

A-" : • - - 2 . - ■ ■ : : \ p 



5(2 -ia)b 
17 — uj = 



(1 - ia)(2 - ia) 

16-30ia-15a 2 ,, 4 -2 + 3ia , 9 

6fc 4 H pk 2 . 

16 4 ^ 



(37) 



2. A stationary pulse fZQ Eq. (119)] 



f A = i4 (-**sechifea:)e*'[ aLo e oosh * a: " wt \ 



2(5 - 4ia - a 2 )b 



Aid 



(1 — ia)(2 — ia) 



(38) 



i7 - w = (1 - ia) 2 [(1 - ia) 2 6fc 4 +pfc 2 ] . 



In both solutions, the number of constraints on (?, r, 7) is two (defined by the vanishing 
of the imaginary part of the relations for k 2 and 17 — u>). 

3 Investigation of the amount of integrability 

3.1 Counting arguments based on singularity analysis 

Because the ODE (J5J is assumed nonintegrable, the number of integration constants which can be 
present in any closed form solution is strictly smaller than the differential order of the ODE. Let us 
first compute precisely this difference, an indicator of the amount of integrability of the ODE. The 
technique to do so is just the Painleve test (see Ref. [Jj for the basic vocabulary of this technique). 
Let us present it on the KS example l|22|l . 

Looking for a local algebraic behaviour near a movable singularity xq (movable means: which 
depends on the initial conditions), 



(39) 



one first obtains the usual balancing conditions (here, between the highest derivative and the 
nonlincarity) 



p - 3 = 2p, p(p - l)(p - 2)vu Q + 



0. 



easily solved as 



p = —3, mq = 120^, 



which yields the Laurent series 



i2o,x- 3 -i5fr- 2 + 15( ;X- b V 

4 x 19. 



I3(4pv -b 2 )b 
32 x 19. 2 



oix 1 ), 



(40) 



(41) 



(42) 



from which two out of the three arbitrary constants are missing. These two constants appear in 
perturbation 



u = + eu 



(i) 



e 2 u^ + . . 



(43) 
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in which the small parameter e is not in the ODE i|22[l . The linearized equation around 



^ + & ^+/^+ u(0) V (1) = ' ( 44 ) 
da;' 3 ax z ax J 

is then of the Fuchsian type near x — xq, with an indicial equation (q = —6 denotes the singularity 
degree of the lhs E of 




\iva X - j - q {vdl+u x v )x j+v (45) 
= v{j - 3)(j - 4)(j - 5) + 120^ = u(j + l)(f - 13j + 60) (46) 

(47) 

The local representation of the general solution, 

u(x ,ec + ,ec_) = 120^x _3 {Taylor(x) 

+e[c + x (13+Jv7T)/2 Taylor( X ) 

+ C _ X ( 13 -^ r )/ 2 Taylor(x)]+0(£ 2 )}, 

in which "Taylor" denotes converging series of x, does depend on three arbitrary constants (a?o, ec + , ec_) 
(the Fuchs index —1 only represents a shift of xq). The dense movable branching arising from the 
two irrational indices characterizes |38| the chaotic behaviour, and the only way to remove it is to 
require ec + — ec_ = 0, i.e. e — 0, thus restricting to a single arbitrary constant the analytic part 
of the solution. 

To summarize, let us introduce two notions. 

The first one is trivial. One calls irrelevant any integration constant which, because of some 
symmetry, is always present in any solution. The KS ODE has one such irrelevant integration 
constant, the origin £o of £, and we will systematically omit to write it. The traveling wave 
reduction of CGL3 has two irrelevant integration constants, the origins of £ and ip, etc. 

The second notion is quite an important property of the equation. We will call unreachable 
any constant of integration which cannot participate to any closed form solution. The KS ODE 
has two unreachable integration constants. 

We will also call general analytic solution the closed form solution which depends on the maxi- 
mal possible number of reachable integration constants, and our goal is precisely to exhibit a closed 
form expression for this general analytic solution, whose local representation is a Laurent series 
like (PI). 

The above notions (irrelevant, unreachable) are attached to an equation, not to a solution. Let 
us similarly introduce two integer numbers, attached to a solution, allowing one to quantify how 
far this solution is distant from the general analytic solution. 

We will call deficiency of a closed form solution the number of reachable integration constants, 
excluding the irrelevant ones, which are missing in this solution. In KS for instance, the elliptic 
solution has a zero deficiency (A is arbitrary), and all the trigonometric solutions have deficiency 
one (A is fixed). 

Let us finally define the codimension of a closed form solution of an equation as the number of 
constraints on the fixed parameters (fixed means: which occur in the definition of the equation). 
Thus, the elliptic and trigonometric solutions of KS have codimension one, and the rational solution 
has codimension two. 



3.2 Evidence for unknown solutions 

Computer simulations as well as real experiments (for a recent review, see [33]) sometimes display 
regular patterns in the (x, t) plane, and some of them are indeed described by some analytic 
solution. For the remaining patterns, the guess is that there should exist analytic expressions, to 
be found, corresponding to these patterns. 

For the KS equation Q), one thus observes a homoclinic solitary wave [23 Fig. 7] (p = /(£), £ = 
x — ct, while all solutions known to date are heteroclinic. 
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For the CGL3 equation (JIJ, it has been predicted J3| the existence of a fourth physically 
interesting solution, which is a codimcnsion-one homoclinic hole solution with an arbitrary velocity 
c. 

Tabic |21 gathers the current state of the known solutions for the various nonintegrable equations 
considered in this chapter. In Ref. [T3|, another counting, based on the various possible topological 
structures, is made for CGL3 and provides the same results. 

Table 2: Integer numbers rating the particular solutions of a nonintegrable equation. The vo- 
cabulary (irrelevant, unreachable, deficiency, codimension) is defined in section f3. II The column 
"Available" indicates the number of relevant, reachable integration constants, this is the algebraic 
sum "Order" - "Irrelevant" - "Unreachable" . The last two columns indicate the properties of the 
solutions in their order of appearance section |3 The best solution would be one with deficiency 
and codimension both equal to zero, with the reduction parameters c and u> arbitrary. 



Equation 


Order 


Irrelevant 


Unreachable 


Available 


Deficiency 


Codimension 


CGL3 


4 


2 = £o,¥>o 


2 





0,0,0 


1,2,2 


CGL5 


4 


2 = £o,<A) 


2 





0,0 


1,1 


KS 


4 


1 = 6 


2 


1 


0,1,1,1,1,1,1,1 


1,1,1,1,1,1,1,2 


SH 


8 


2 = £o,¥>o 


6 





0,0 


2,2 



4 Selection of possibly single valued dependent variables 

Whatever be the class of methods to be applied, a prerequisite is to determine a variable whose 
dominant behaviour is single valued and which satisfies some algebraic ODE (or more generally 
PDE). This is the case for KS since the solution p of i|40|) is integer, but in CGL3, CGL5, or SH 
this is the case of neither (A, A), nor argyl, and for CGL5 this is not even the case of \A\. Indeed, 
considering CGL3 for instance, the dominant terms are 

pA xx + q\A\ 2 A, (48) 

and one easily checks that \A\ generically behaves like simple poles. Let us therefore define the 
dominant behaviour of the two fields (A, A) as 

A ~ A oX ~ 1+ia , A ~ l^- 1 -"*, A £C, aeK, (49) 

in which (Ao,a) are constants to be determined. The resulting complex equation (equivalent to 
two real equations) 

(-l + ia)(-2 + ia)p+\A Q \ 2 q = 0, (50) 

is precisely the one artificially introduced earlier as (|llf> and solved in 115|) . the convention that 
Aq is real being allowed by the phase invariance of CGL3. The same applies to CGL5 

A ~ A oX - 1/2+ia , A ~ A oX - 1/2 - ta , (51) 

see (SH and and to SH, 

A~A oX ~ 1+ia , A-Aox- 1 - 1 *, (52) 

see PS|. 

In all three examples CGL3, CGL5, SH, the variable M = \A\ 2 satisfies an algebraic ODE, which 
can be built by elimination of argA, and it has a single valued dominant behaviour (respectively 
movable double poles, simple poles, double poles). Moreover, again for CGL3, CGL5, SH, for all 
the solitary wave solutions which are known to date (see section [2J, this variable M is represented 
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by quite simple mathematical expressions, namely either polynomials in one elementary variable 
t, Eq. (|16fl . which satisfies a Riccati equation 

— t(z) = 1 — r 2 , t = tanhz, (53) 
dz 

or (source solution of CGL5, Eq. I|34|ll polynomials in two elementary variables (a, r) which satisfy 
a projective Riccati system jS] 

— t = 1 — r 2 — Uqct, — a = —<jt, a 2 — t 2 — 2/iocr + 1 = 0, (54) 
dz dz 

in which /Uq is a constant, and whose solution can be expressed as 

sinh z sinh ka 



-, IM) = cothfca. (55) 
cosn z t cosii ku cusn z -|- cusii /tu 

When /io(Mo — 1) = 0, the class of polynomials in (a, r) degenerates to the class of polynomials in 
(sech, tanh) . 

Therefore M = \A\ 2 will be our best choice to search for closed form solutions of CGL3, CGL5, 

SH. 

Remark. Despite the multivalued dominant behaviour of the complex amplitude A of CGL3 
and SH, one can define two variables with a single valued dominant behaviour. In this complex 
modulus representation [§| 

A = A Z(0e i ^- ujt { A = A Z(Oe- i ^- ajt \ (56) 

with Z complex and $ real, the dominant behaviour is 

Z~ X -\ Z~x~\ &~a X -\ (57) 

and the truncation of (Z, Z, $') might prove to be much more economical than that of M. All the 
solutions listed in section [21 for CGL3, CGL5, SH have been written in this representation. 



5 On the price to obtain closed form expressions 

Let us now give some details on the distinction between the two main classes of methods outlined 
in the introduction. 

In the first class of methods, one gives as an input some class of expressions /(£) (for instance 
polynomials in sechfc£ and tanhfc£), and by a direct computation one checks whether there indeed 
are some solutions in the given class. We will call for shortness these methods sufficient, because 
they for sure miss any solution outside the given class, e.g. for the ODE 

12M 2 - - M' + 36M 4 M 2 + - = 0, (58) 

2 7 2 2 

its solution since it is rational in tanh fc£, 

tanh(£ - Co) 
2 + tanh 2 (£-£o) 

In the second class of methods, the search for first order autonomous subequations J7J) 
requires no a priori assumption at all, and, from the classical results recalled in Appendix, the 
knowledge of the first order subequation is indeed equivalent to the knowledge of the explicit 
expression |JBJ|. As opposed to the previous methods, which are "sufficient" as said above, the 
proposed method can be qualified as "necessary" . 

The difference between the two classes of methods is obvious: the class of expressions /(£) is an 
output of the second method, while it is an input of the first one. This is why the second method 
can find, if they exist, not only some but all the solutions which are elliptic or trigonometric. 

Remark. The cost of the method of first order autonomous subequations is an increasing 
function of the positive integer m occuring in H168(l . but to, which is an input of the method, is 
not bounded. Indeed, any rational function u = Pjv(tanh(£ — £o))/-Fi>(tanh(£ — £o)) satisfies an 
ODE (7} of order one and degree max(iV, D). By considering only some differential consequence 
of this ODE, one cannot guess the correct value of m in advance. 



M = ™ . (59) 
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6 First class of methods: truncations 



After having selected, as indicated in section 0] dependent variables with a single valued leading 
behaviour, the methods called truncations consist in defining for each such dependent variable 
some single valued closed form class of expressions, then in checking whether there exist solutions 
in that class. 

The class of expressions to choose as an input depends on the number of families of movable 
singularities of the considered dependent variable. Thus, the field u of KS has only one family, 
i.e. one value of uq, while the field M = \A\ 2 of CGL3, CGL5 or SH has respectively two, four, 
and four families. Let us start with the simplest class of expressions. 

6.1 Polynomials in tanh (one-family truncation) 

The class of polynomials in tanh(fc/2)£ is the most frequently encountered class of closed form 
solutions of autonomous PDEs. This fact is the direct consequence of a quite remarkable property. 
Indeed, by a result of Painleve, the variable r in l|53|) is the unique variable to be at the same time 
single valued and closed by differentiation: if u is such a polynomial, 

« = E u ^ 3+p ' = \ tanh \z> £ = x - ct > ( 6 °) 

3=0 

the lhs E(u, x, t) of the equation of the PDE is also such a polynomial, 

E = Y,E jX j+ \ (61) 

3=0 

and its identification to the null polynomial 

Vj :E j = 0, (62) 

generates the smallest possible number of determining equations Ej = 0. As compared to the 
Laurent series (|42|l . the series (|60|l terminates, hence its name of truncation. 

The truncation l|60|l involves only one value of uq, it is called for this reason a one-family 
truncation. Let us give a few examples. 

6.1.1 One-family truncation of the KS equation 

The symbols uq and p denoting the leading behaviour of the ODE (1221) , the truncation l|6(J|) defines 



— q + 1 = 7 determining equations, the first four being 

u 2 

E Q = -60isu + y = 0, (63) 

Ei = 126m + (no ~ 24j/)ui = 0, (64) 

57 1 
E2 = — 3/iUo + -^k 2 vuo + 6bui + — u x + (uq — &v)u2 = 0, (65) 

g 

E% = —-bk 2 uo — 2/iUi + 10k 2 vui + 2bu2 + U1U2 + u^u^ = 0. (66) 



The structure of this kind of algebraic determining equations is always the same: one algebraic 
equation for uq (j = 0), followed by — p equations linear in Uj,j = 1, • • • , —p. Equation j = has 
already been solved, see ()41[) . and the next equations j = 1, . . . , — p have the same solution Uj as 
in the infinite Laurent series 142(1 . The truncated expansion (|60l) then evaluates to 

u = T> Log rp + constant, (67) 

in which T> is the singular part operator defined in l(27|l from the Laurent series, and ip is the 
logarithmic primitive of x 1 an entire function defined by 

r - ^ = 0, (68) 
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whose value can be chosen without loss of generality as 

V> = cosh-f. (69) 

After the operator V has been computed, the two equations l|67(l . (|fiH|> are an equivalent way 
of defining a one- family truncation, much more elegant than with (|60|) . 1)53(1 . 

The remaining — q + p equations j = —p + 1, . . . , — q are algebraic in k 2 and the parameters 
appearing in the definition of the equation 1221) (one says the fixed parameters), 

E A = --b 2 k 2 + —ii 2 + —b A u~ 2 - —b 2 ^- 1 + 40fcV - 76JfeV = 0, (70) 
2 19 P 304 38 P P v ' 

E 5 = b(b 2 - 16,H ($k 2 + ^b 2 v- 2 - = 0, (71) 

E 6 = 32A + 3vu k e + 4(6«i - z/u 2 )fc 4 + 8fc 2 /w 2 + Wuj = 0, (72) 
and they admit only the six solutions listed in Tabled 

6.1.2 One-family truncation of the real modulus of CGL3 

Whatever be the chosen representation (couple (M, <p) 1 (Z,Z,Q), etc), the CGL3 equation has 
more than one family, sec (|15fl . therefore any one- family truncation only captures part of the 
whole singularity structure and cannot yield the general analytic solution. Nevertheless, as already 
noticed in section l2~D Ea. ((16J) . the one- family truncation of M = \A\ 2 must provide at least the 
three currently known solutions. Let us perform it. 

The field M has two families of movable double pole-like singularities 

M-^|^X- 2 (l + fx + 0(^)), (73) 

with singular part operators T>± equal to 



In its elegant definition, the one-family truncation, 



M = T>± Log -0 + m, 

S k 2 (75) 

i/)" + —if) = 0, S — — — — constant, 



transforms (|19|) into the truncated Laurent series 

14 

5>^- 14 = 0, (76) 

and one must solve the 15 real determining equations Ej = for the two constant unknowns 
(S, m) and the five parameters d r , di, g r , gi, csi occurring in l|19f) . By construction of T>±, equations 
Ej = 0, j = 0, 1, are identically zero. 

To avoid carrying heavy expressions, let us make the following nonrestrictive simplification. 
Out of the five parameters d r ,di,g r ,gi,csi of the ODE (|19fl . only three are essential (g r ,gi,c, 
equivalent to 7, u>, c). Indeed, p and q (i.e. d r + idi and s r — isi) can be rescaled to convenient 
numerical values, such as 

p = — 1 — 3i, q = 4 — 3i, 
, 1 , 3 1 3.9 

d r = — , a, = — , s r — , s, ■ — , A = — . 77 

2' 2' 10' 10' 2 y ' 
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Choosing the + sign in J7SJ), one has T> + = 4(— d 2 — (c/10)d x ). As seen from the first few 
determining equations, 



16g r = 0, 



57 

E 2 = —c 2 + 156c 2 + 13k 2 + 4 
E 3 = (-'^c 2 - 432c 2 - 28k 2 - 10,,, -Is,/, j« = 0. 



(78) 
(79) 



the resolution presents no difficulty. In particular, after solving the equations numbered j = 
0, . . . , 6, all the remaining equations are identically zero, a fact which indicates a high redundancy 
in these determining equations, which are therefore not at all optimal. In the CGL3 case properly 
said lm(p/q) ^ 0, for each sign in l|74|) one obtains three solutions, 



M = -2 



T — 

2 4 



c \ 2 / c y 
2o) + (lo) 



2 



c c 
T_ 20 ~ 5M 



/ C \ z 4 



M 



-2 r 



9 



3c^ 
50 

T, 



0. 



fc = 2g r , c = 0, #i = 0, 

2 



M 



2 T± 



k 



2 



20' 



,9r = 0, gi - — = 0, 



and 



= 4 



c 

20 

2 2 



/ c \ z I 

VloJ + 3 3r ' 5l ~ 9r + 



3V 
80 



c£r _ _ _c c_ 

2 ~~ T 20 5M 

= 0, 



T 2 - 



2 2 

k " = g5r, c = 0, 3g; - g r = 0, 
2 



M = 4 r ± 



2r 



10' 



10 



fc2 = 77T . 5r = 0, gi - — = 



50 



These solutions are identical to those listed, in the same order, in section |2~T1 

6.1.3 One-family truncation of CGL3 in the complex modulus representation 

As already outlined at the end of section^] the one-family truncation of (Z, Z, $') 

Z = x' 1 + X + iY, 
Z = x 1 + X — iY, 
$ = aLogi/' + K£, 



with the gradient definitions 



(Log^)'-X- 1 , 
k 2 

x' = i - T x 2 , 



puts the lhs of Eq. in the form 



(80) 



(81) 



(82) 



(83) 



(84) 



(85) 



(86) 



(87) 



3-3 _ 



0. 
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thus generating four complex determining equations Ej = 0, (i.e. eight real, to be compared with 
the fifteen of section l6\l.2[) . These equations must first be solved as a linear system on C, as follows 
HI Appendix A]. The first equation E = 0, identical to l(TT|) . is linear in p and q, let us solve it 
for q, 

q = -(1 -ia)(2-ia)A Q 2 p. (89) 
The next equation j — 1 is then linear in K, X, Y, c, let us solve it for instance for K, 

K = (Si + a)X -Y + (90) 
2p 



The equation j = 2, linear in 7, ui, k 2 , is solved for (17 — oj)/p 

2 

ij — U) 
P 



= [ Yp ) + - (1 - - (2 - 3t 



(91) 



and the advantage of this pivoting elimination is that the last equation j '• — 3, which does not depend 
on (g, if , 7) by construction, is also independent of (p, c, a;, j4o). It only depends on (X, Y, a, k 2 ), 
and it factorizes as 



£3 = (2X - aY)(4{X + iYf - k 2 



0, 



(92) 



thus defining two solutions on C. 

Finally, considering now the system H89|) - (|92|l for the real unknowns or parameters (Aq, a, K, c, X, Y, 7, u>, k 2 ), 
it is quite easy to obtain the three solutions listed in section |2~T1 

6.2 Polynomials in tanh and sech (two-family truncation) 

The class of polynomials in tanh and sech 



'-p-i 



a,j tanh kt; + I fry tanh fc£ sechfc£, 



(93) 



3=0 



can equivalently be represented by the class of powers of tanh ranging from p to — p |31| , 



-2p 
3=0 



k k 

+ X" 1 = g tanh 2^' ^ = X ~~ Ct ' U U -2p ^ °) 



because of the elementary identities 9 
1 



tanh z — 



tanh z 



-2i sech 



7f 

2z + ,- 



tanh x 



1 



tanh z 



2 tanh 



2* + .- 



(94) 



(95) 



A solution in this class can only exist for ODEs admitting at least two families with the 
same p. Indeed, if for this p there exists only one value of Uo, only the second combination 
tanh +1/ tanh = 2 tanh can contribute. For instance, the KS equation cannot admit such a solution. 

More generally, the class of polynomials in r and a defined in l|55|) , 



V 3=0 

is equivalently defined as |9j Appendix A] 



-p-i 



J2 bjT* ] a, (o_p,6_ p _ 1 ) 9 fe(0,0), 

3=0 



(96) 



u = T>\ Log ijii + T>2 Log ip2 + Tn, m = const, 

S S k 2 

4>" + g^ 1 = °' ^2 + 2^ 2 = °' S = ~T = constant ' 



s 



-01 "02 



(97) 
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In this writing, which is the natural extension of (|75|l to two families, the linear operators T>i and 
T>2 are the singular part operators of two different families, the entire functions and tp2 obey 
the same second order linear equation, but with a different choice of the integration constants, 

k k 
ipi = cosh — (£ + a), ip2 — cosh — (£ — a) , /iq = cothfca. (98) 

The case /xo(/Xq — 1) = reduces to the class of polynomials in tanh and sech. 
The practical implementation is the following. 

1. For the class of polynomials in tanh and sech, one puts the lhs E(u) of the nonlinear ODE 
under the same form as u, 





-2p 




u = 




-2p + 0, 




0=0 








k 2 


x' = 


= i + fx 2 ,s = - 


2 ' 




-2q 




E = 








0=0 




Vj: 


Ej = 0. 





(99) 



and one solves the set of — 2q + 1 determining equations Ej = 0. 



(100) 



2. For the class of polynomials in r and a defined in l|55|l . under the assumption H97J1 . the 
lhs E(u) is first expressed as a polynomial of the two variables ip'j/t/jj,j = 1,2 



-9 -q-k 



fe=0 1=0 



V>2 



(101) 



which further reduces, thanks to the third line of (|97() . to the sum of two polynomials of one 
variable, 




(i) 



(102) 



One then requires the vanishing of the — 2q + 1 determining equations 



£0 = 0, E 



(i) 



0, E. 



(2) 



0, j = 1, 



(103) 



As an example, let us apply this to the ODE 

dtt x ~ 



E(u)= (^) -a 2 (« 2 -6 2 ) 2 +c = 0. 



(104) 



It admits two families, with singular part operators T>\ = a. d^,T>2 = —a. 1 . The relation 

T> 2 = -Pi implies Ef ] + Ef ] = 0, j = 1,2,3,4, and only 5 out of the 9 determining 

equations are linearly independent. Moreover, by construction of the singular part operators, the 
two equations j ' = 4 are identically satisfied. The next equation j ' = 3 



(i) 



-2a fc/io — 4a m = 0, 



(105) 



is solved for m. Then, the equation j = 2 



E. 



(i) _ 



26 2 + a~ 2 k 2 - -a- 2 (kfi ) 2 = 0, 



(106) 
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is solved for k 2 , considering kfi as a single variable. The remaining system 

E[ l) = (~2b 2 + i a - 2 (fc Mo ) 2 ^ (fc M0 ) - 0, (107) 

E = -a?b* + c + \b 2 {k^) 2 - l a - 2 (fc Mo ) 4 = 0, (108) 
2 lb 

admits two solutions. The first one c = 0, (fc/io) 2 = (2ab) 2 corresponds to a factorization of the 
equation E(u) — into two Riccati equations and therefore must be rejected. The second one 

fc/io = 0, a 2 b 4 - c = 0, (109) 

defines a solution, provided the indicated constraint on the fixed parameters (a, b, c) is satisfied. 
This solution 

Mo = 0, m = 0, fc 2 = -2 M ) 2 , - = «-^Log^§|±^ (110) 

d£ cosh(£;/2)(£ — a) 

is nothing else than u — i(k/a) sechfc£, using the relation (1q = cothfca. 

Indeed, as opposed to the function tanh, which satisfies an ODE admitting only one family of 
movable singularities (the Riccati equation), the function sech (or more generally its homographic 
transform a) satisfies a first order second degree ODE 

sech' 2 + sech 4 - sech 2 = 0, (111) 
which admits two families of movable simple poles with opposite residues 

sech^-^-i^-Cor 1 . (112) 

6.2.1 Two-family truncation of the real modulus of CGL3 

M admits exactly two families, so its two-family truncation is quite appropriate. The two singular 
part operators are, for each family, defined in l|74l) . with 2?i = T> + ,T>2 = T>-. The assumption 
(|97|l . with p = —2, q = —14, transforms ljl9|l into the sum Hl()2(l of two polynomials of one variable, 
and the four equations j = 14, 13 are identically zero, by definition of T>±. 

For p, q, 7 arbitrary, the resolution of the 25 remaining determining equation is impossible to 
carry out by hand. But the hand computation becomes possible with the generic numerical values 
(|77jl . First, the system j = 12 

E& = — c 2 - 780m - 520fc 2 + 78cfc^ + 390(fc/t ) 2 + 80. 9l + 320. 9r = 0, (113) 
5 

q 

E^ = -c 2 + 120m - 40fc 2 - 24ckfi + 120(kfi ) 2 + 20g t - A0g r = 0, (114) 



is solved as a linear system for m and k 2 . The next system j = 11 



= _^ c 3_^ c 2 fe/io + ^£ c(A . Mo) 2_ 312(fcMo) 3 



39 , 254 2 , 468 

C C kiln H 

125 25 M 5 
K~ c 9i ~ 168/c/i 5i 5~ C5r ~ 48fc Moffr = 0, (115) 

E{f = -^-c^ + ^^k^-Z^k^f + lb&ik^f 
87 

— r c 9i + 84fc/i ffi - 2cg r + 2Ak^ a g r = 0, (116) 
5 

is linear in (g r ,gi), with a Jacobian J ~ c(3c — 5fc//o)- For the first subcase J ^ 0, after solving 
for (g r ,gi) as functions of (c,k/in), the next system j = 10 only depends on k^jc and it admits 
no solution. The discussion of the second subcase J = leads to the conclusion, only using the 
next system j = 10, that no solution exists. An identical result is achieved for arbitrary values of 
(p, q, 7) using computer algebra. 
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This unfortunate situation is exceptional, and only reflects the difficulty of CGL3. Should such 
a solution exist, it would have the form 



M 



— > 11 9<i r 2 
— ^-tanh+ci sech+ — ^ tanh +C3tanh+C4, 



2d 2 



2d? 



(117) 



and the constraint C3 = could define a homoclinic hole solution, just like the (yet analytically 
unknown) one of van Hecke |13) . 

6.2.2 Two-family truncation of CGL3 in the complex modulus representation 

Let us denote (A ,a) and (^2,012) two different solutions of l|15|l. 

In the complex modulus representation (|56J) . the two-family truncation of (Z, Z, $) is defined 
as Appendix A], 



A = (A (8 S Log MO + X + 1Y) + A 2 8^ Log MO) 



3 <[-wt+*(0] j 

A = (A (d 6 Log MO +X-iY) + A 2 ds Log MO) e~^ t+ ^, 
$ = a Log -01 + a 2 Log i\) 2 + Kt;, 



(118) 



with the definitions for the derivatives of (ipi,ip2) given by the last two lines of l|97|l . The lhs of 
Eq. then takes the form (|l(J2[l with q = —3, and one solves the seven complex determining 
equations as a linear system on C, similarly to what has been done in section lbM.3l 
From the two equations j '• = 3, 



E ( 3 1} = A ((1 - ia)(2 - ia)p + A 2 q) = 0, 



E. 



(2) 



A 2 ((1 - ia 2 ){2 - ia 2 )p + A 2 2 q) = 0, 



and the two relations implied by (TSj, 



4 12 di . 2 

a = -^A , a 2 = —A 2 . 



(119) 
(120) 



(121) 



one first proves that the only possibility is A 2 = —Aq, therefore the two complex equations j ' = 3 
are solved as 



A 2 = —A , a 2 = a, q — -(1 - ia)(2 - ia)A a 2 p. 
At the level j = 2, the symmetric combination 

E [ 2 1] + E { 2 2) = p(l - ia) 

is solved for the two pieces of information 

X = 

then the antisymmetric combination 



(ia — 3)X — iY + ( ia — — ) k/io 



~ kfi , Y = ^akfio, 



= 0, 



- E {2) = p(l - ia) 



2p 



- K 



= 0, 



is solved as 



K 



2p 



(122) 



(123) 



(124) 



(125) 



(126) 



At the level j = 1, the symmetric combination is identically zero, and the antisymmetric combi- 
nation is solved for (ij — uj)/p (we omit the expression). The remaining equation 



E = Mo [(2 + (a 2 - 2)fi 2 ) + ia{-4 + (a 2 + 4) M g)] = 



(127) 
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admits as only solution /i = 0. Therefore, one obtains the unique solution 

A 2 = -A , a 2 = a, q= -(1 - ia)(2 - ia)A~ 2 p, (128) 

X = 0, Y = 0, A t = 0, K = (129) 

2p 



, +{l-iaYk\ (130) 
P \2pJ 

The value of K implies cpi = 0, and the case c = represents the homoclinic pulse (|13|) . 
6.3 Polynomials in p and jp' 

A class of polynomial elliptic functions can be defined for instance with the Weierstrass function 
p and its derivative ^| , 

K-i>-3)/2] \ 

E MO' P'K). (131) 
i=o / 

Since p admits only one family, such a solution may exist for any ODE. It will be quite useful to 
take advantage of the value of the singular part operator of p(£), 

]2 




V = -— 2 - (132) 



For KS, the assumption to seek for solutions in the above class 

u = c Q p' + cip + c 2 , c ^ 0. (133) 

together with the knowledge of the singular part operators (|27|l and (|132fl . first yields the correct 
values of cq and c\, 

u= -60vp' -15bp + c 2 , (134) 

a truncation which defines the four determining equations 

b 2 — 16/ii/ = 0, 
fyu + 4vc 2 = 0, 

bc 2 - 720v 2 g 2 = 0, ( 135 ) 
A +^c 2 2 + ^- b 2 g 2 + 30nvg 2 - 1080^ 2 .g 3 = 0. 

Their unique solution is (|23|1 . 

For CGL3, the assumption that M be in this class 

M = a 2 p + c 2 , a 2 t^O, (136) 

generates 10 determining equations, in the parameters and unknowns (a 2 , c 2 ,g 2 , g^ d r , di, s r , csi,g r , gi). 
The equation with the highest singularity degree 

(Ka 2 ) 2 - 9a 2 d r - 18) (2 + a 2 d r ) = (137) 

in which the vanishing of the second factor is forbidden, is first solved as a linear equation for d r 

dr = w^. (138) 

9a 2 

The next equation yields csi = 0. The next equations are successively solved for g r ,g2,g3, and 
the elliptic discriminant g 2 — 21 g 2 is then divisible by the unique remaining determining equation. 
Therefore, one finds as unique solution the pulse (|13|l . 

Finally, let us mention the Ansatz made for CGL5 |25l 13"] 

A = a(x)e^- 2aLo & a ^' ujt \ (to, a, a) eU, (139) 

which sets an a priori constraint between the amplitude and the phase (similar to that made for 
CGL3 in Ref. [5]), together with the assumption that a 2 obeys a first order second degree elliptic 
equation. This allows one to retrieve (|3"4l in the particular case r = 0, c = 0. 
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7 Second class of methods: first order subequation 

The requirement that the solution ||JjJ be shared by the iVth order ODE (JjJJ and the first order 
ODE Q characterizes, as recalled in the Appendix, the singlevalued expressions / as being elliptic 
or degenerate elliptic (i.e. trigonometric or rational), i.e. the class 

u = R(p',p) ^i?(e fc «) — >fl(0, (140) 

in which R denotes rational functions and — ► denotes the degeneracy. This class contains all the 
classes considered in previous sections (polynomials in tanh, in (tanh, sech), in (a, r), in (p, p')), 
but it also contains in addition expressions like l(59|l . 

The algorithm to obtain all the elliptic solutions combines two pieces of information: 

1. a local one, a Laurent series representing the largest analytic solution of the N-th order ODE 
near a movable pole-like singularity 

2. a global one, the necessary form (|168|l of Q. 

by requiring that the Laurent series satisfies the first order subequation l|168fl . 

This provides the explicit form of the first order subequation F(u, v!) — 0. Then one computes 
the solution u — /(£ — £o) from this equation F = 0. 

The successive steps are [23 Section 5] . 

1. Choose a positive integer to and define the Briot and Bouquet first order ODE 

m [(to— k)(p-l)/p] 

F(u,u')=J2 E a jtk u j u' k = 0, a , m ± 0, (141) 

fc=0 j=0 

in which [z] denotes the integer part function. The upper bound on j implements the con- 
dition m(p — 1) < jp + k(p — 1) that no term can be more singular than u' m , identically 
satisfied if p = — 1. The polynomial F contains at most (to + l) 2 unknown constants 

2. Compute J terms of the Laurent series, with J slightly greater than the number of unknown 
constants a^fc. 

u = * P ( E + °(X J+1 ) ] , X = f - Co, (142) 

\3=0 

where p is —3 for the KS equation (|22H . —2 for the variable |A| 2 of CGL3, etc. 

3. Require the Laurent series to satisfy the Briot and Bouquet ODE, i.e. require the identical 
vanishing of the Laurent series for the lhs F(U, U') up to the order J 

F = x D (j2 + °(X J+1 ) j , D - m(p - 1), (143) 
Vj : Fj = 0. (144) 

If it has no solution for Oj,fc, increase m and return to first step. 

4. For every solution, integrate the first order autonomous ODE 1)141(1 . 

Let us give two examples. 
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7.1 First order autonomous subequations of KS 

The Laurent series of is 

In the second step, the smallest integer m which allows a movable triple pole (p = —3) in (|141fl 
is vri — 3. With the normalization 0,0.3 = L the subequation contains ten coefficients, which are 
first determined by the Cramer system of ten equations Fj = 0, j = : 6, 8, 9, 12. The first few are 

F = -9a , 3 + 4fVa 4 , = 0, (145) 

Fi = 96a , 3 + 12i/ai, 2 - 806z/a 4j0 = 0, (146) 

F 2 = (21206 2 + 25 60^)^4,0 - (1056 2 + 144/^)a , 3 - 5326^ai, 2 - 608^ 2 a 2 ,i = 0, (147) 
F 3 = (5b 2 + 72(j,v)ba , 3 + (1376 2 + 240fiv)vai, 2 - (44 26 2 + 2656/^)6m 4 ,o 

+6086i/ 2 a 2 ,i + 608z/ 3 a 3 ,o = 0. (148) 

The remaining infinitely overdetermined nonlinear system for (v, b, (i, A) contains as greatest com- 
mon divisor (gcd) 6 2 — 16fiv (see Eq. (|23f) 'l. which defines a first solution 

b 2 1 / . b \ 2 f , b \ 9 / 2 156 6 10A\ 2 n 
= 16, (u'+— Us) \u' -—u s ) +—[u 2 s + T ^j- I + —J =0, 

(149) 



(iv V 2v J \ Av J 4Qv \ s 1024v 

3b 3 



32v 2 



After division by this gcd, the remaining system for (y, 6, (l, A) admits exactly four solutions (stop- 
ping the series at j = 16 is enough to obtain the result), namely the first three lines of Tabic ^ 
each solution defining the same kind of subequation, 



b = 0, 

. 180/i 2 \ 2 / , 360/i 2 \ 9 / 2 30it . 30V\ 2 , . 

b = °> u + wA u+ ^ u+ wj =0 ' (151) 

52 144 563 f / fo V 9 4 n ,, Bn s 

-„ = ^ u ° = u -u^> \ U+ T V U °) + 40^ = ' (152) 

b 2 256 456 3 



/if 73 ' 2048^ 2 ' 

b V ( , 6 \ 9 / 2 5o 3 56 2 , . 

-I u s it H u s H u{ -\ ~u s H u = 0, 

8v S J \ 2v s 40i/ V 1024t/ 2 128^ ' 



(153) 



In order to integrate the two sets of subequations (jl49[l . H150[) - H153|) . one must first compute 
their genus , which is one for (|149|l . and zero for l|150(l - H153|l . Therefore (|149|) has an elliptic 
general solution, listed above as l|23l) . and initially found |12M17| by other methods. 

As to the general solution of the four others i|150[l - ill53|) . this is the third degree polynomial 
(EJ in tanhf(£-£ ). 

These four solutions, obtained for the minimal choice of the subequation degree m, constitute 
all the analytic results currently known on Ij22(l . For m — 4, no additional solution is obtained |41| . 
The computation for m = 5 is in progress. 

7.2 First order autonomous subequations of CGL3 

We consider the variable M = \A\ 2 , i.e. p = —2. The smallest value of m is then 2. With the 
numerical values 1|77[) . the two Laurent series are 

M- =X- 2 (-2+ C -X+(^-^-^)x 2 + 0( X 3 )) , (154) 



5^ V 3 6 200 / 



^^For instance with the Maple command genus of the package algcurves 1151 . which implements an algorithm of 
Poincare. 
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The existence of two Laurent series, rather than just one, is a feature which the subequation 
must also possess, and this has the effect of setting the lower bound to m — 4 instead of 2. Indeed, 
the lowest degree subequations 

F 2 = M' 2 + M\a hl M + a ,i) + a 3j0 M 3 + a 2 , M 2 + a h0 M + a ,o = 0, (156) 
F 3 = M' 3 + M' 2 {a ia M + a , 2 ) + M'(a 3 ,iM 3 + a 2A M 2 + a hl M + a ,i) 

+a 4 , M 4 + a 3i0 A/ 3 + a 2 , M 2 + a lj0 M + a 0:0 = 0, (157) 

have the respective dominant terms M' 2 + a^ oM 3 and M' + a^ iM' M 3 , which define only one 
family of movable double poles. 

Let us nevertheless start with m = 2, for which (|156|l can only be satisfied by one series, 
e.g. (|154J) . thus preventing the full desired result to be obtained. The six coefficients aj.k of (|156|) are 
first computed as the unique solution of the linear system of six equations Fj = 0, j = 0, 1, 2, 3, 4, 6. 
Then the J + 1 — 6 remaining equations Fj = 0,j = 5, 7 : J, which only depend on the fixed 
parameters (g r ,gi,c), have the greatest common divisor (gcd) 3<?i + 2g r + 3c 2 /50, and this factor 
defines the first solution f l|158(l below). After division par this gcd, the system of three equations 
Fj =0,j = 5, 7, 8, provides two and only two other solutions, see Ijl59|l and (|160|) below, with the 
respective constraints (c = 0, gi = 0) and (g r — 0, 50^ — c 2 = 0), and all the remaining equations 
Fj =Q,j > 9, are identically satisfied. 

Therefore, with this lower bound m — 2, one already recovers all the known first order sube- 
quations. These are, with the series (I154|l . 

c c 3 \ 2 / c 2 \ / c 2 2 \ 2 3c 2 

M ' + 5 M+ 25o) +2 { M+ V0){ M ~V0~-^) =0,3 ft + 2* + - = 0,(158) 

M' 2 + 2(M - g r )M 2 = 0, c = 0, 9l = 0, (159) 
(M' + |Af) + 2A/ 3 = 0, g r = 0, 9i - ^ = 0. (160) 

Finally, for each of the three subequations, the fourth step finds a zero value for the genus and 
returns the general solution as a rational function of e a ^~^°\ which basic trigonometric identities 
then allow to convert to the second degree polynomials in (fc/2) tanhfc(£ — £o)/2 listed in H80f) ~ H82[) . 

Similarly, with the other series Q155[) . one obtains 

c c 3 \ 2 / c 2 \ / c 2 2 \ 2 3c 2 

M'+-M ( M )[M+ g r ) =0, 3oi-o r + — = 0,(161) 

5 500 / V 100 M 100 3 y / ' y y 80 ' v ; 



M' 2 - M lM - ~ 9r ) = 0, c = 0, 3 3 , - g r = 0, (162) 

Af + 1m) 2 - A/ 3 = 0, g r = 0,gi-SL = 0. (163) 
5 / 50 

With the correct two-family lower bound m — 4, which corresponds to 18 unknowns aj.k and at 
least 24 terms in the series, we have checked that there is no solution other than the above three. 
This situation is quite similar to the absence of solution in the class i|117|l . and it just reflects the 
difficulty of the CGL3 equation. 

The case m = 8 (60 unknowns aj^k and at least 66 terms in the series) is currently under 
investigation but preliminary results seem to indicate the absence of any new solution, and we are 
now automating the computer algebra program in order to handle much larger values of m. 
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7.3 Domain of applicability of the method 

As we have seen, the subequation method contains the truncation methods and its cost is minimal 
since the main step is a linear computation. 

The two key assumptions behind this "subequation method" are, 
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1. a Laurent series should exist, 

2. a first order autonomous algebraic subequation should exist. 

Its best applicability is therefore nonintegrable 7V-th order autonomous nonlinear ODEs admit- 
ting a Laurent series which only depends on one movable constant, such as the CGL3 ODE (|19|l 
or the traveling wave reduction ill' 21) of the Kuramoto-Sivashinsky equation [81 141j. 

Two examples of inapplicability are 

1. the Lorenz model, in which the Laurent series generically does not exist and has to be replaced 
by a psi-series |36| . 

2. the autonomous ODE u'" — \2uv! — 1 = 0, which admits the first Painleve transcendent as 
its general solution, a case in which no first order subequation exists. 



8 Conclusion 



How do these two classes of methods (truncations, first order subequations) really compare, inde- 
pendently of the amount of computation involved? 
Let us first recall a preliminary, classical result. 

The class ()60|) of polynomials of degree — p in tanh obeys a first order equation of degree 
m = —p. For instance, given the polynomial 



u = tanh 2 +2a tanh +b, 
this amounts to eliminate tanh between the two algebraic equations 



tanh +2atanh+fo — u 
2(tanh+a)(l - tanh 2 ) • 



= 0. 

u' 



0. 



(164) 



(165) 



which results in 2 
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2 2a 


-u' 
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2a -u' 







= (u' 


- Aa(u - 


-b + 


a 2 )) 2 -A{u- 


- b + 2a 2 



l) 2 (u- b + a 2 ) = 0, 



(166) 
(167) 



an equation with degree m = 2 = — p, having genus zero. 

Similarly, the class of polynomials of global degree — p in (tanh, sech) or (er, r) obeys a first 
order ODE with degree m — —2p. 

Last, the class (fl31l> of polynomials of (p, p') of singularity degree at most equal to p obeys a 
first order ODE with degree m = —p. 

Therefore, given a value of p (the singularity degree of the ODE) and a truncation considered 
in sectional there exists a value of m (either — p or — 2p) at which the result of the truncation can 
be found by the method of first order subequations. 

Conversely, given a value of m (the degree of a first order subequation), the class of solutions 
of the method of first order subequations is made of the rational functions (of (p, p') or of e k ^ i.e. 
(fc/2) tanhfc£/2), a class richer than the polynomials. 

This proves the identity of the two classes of methods, provided the truncations assume rational 
functions instead of polynomials. 

However, from the practical point of view of the amount of computation involved, the increasing 
order of difficulty seems to be 

1. Truncations of polynomials. 



2 This formula, due to Sylvester, expresses the resultant of two polynomials of degrees m and n as a determinant 
of order m + n. 
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2. First order subequations. 

3. Truncations of rational functions. 

How does this compare with the approach of Chow (see e.g. [5]) to find solutions of PDEs in 
terms of elliptic functions? To be definite, let us consider a PDE in (x. t). The solutions which do 
depend on both x and t (i.e. which do not satisfy some ODE) are richer than those here described. 
As to the solutions of the solitary wave type f(x — ct), the method of Chow belongs to the first 
class of methods, i.e. it may or it may not find the most general elliptic solution which exists. 
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9 Appendix. Classical results on first order autonomous 
equations 

The following results were mainly obtained by Briot and Bouquet, Fuchs, Poincare and put in final 
form by Painleve (22 pages 58-59]. 

Theorem. Given an algebraic first order autonomous ODE Q, whose general solution is 
therefore JBJ, the following properties are equivalent. 

1. Its general solution is singlevalued. 

2. Its general solution is an elliptic function, possibly degenerate. 

3. The genus of the algebraic curve is one or zero. 

4. There is equivalence between the knowledge of / and that of F, 

5. There exist a positive integ er ?7i cLiid (jn H- 1) complex constants flj,fc) with. fto,m 

7^ 0, such 

that the polynomial F of two variables has the necessary form 

m Ira— 2k 

F(u, u') = J2 11 !■>■" ' "' L = 0, a , m ? 0. (168) 

k=0 j=0 

6. If the genus is one, there exist two rational functions Ri,R.2, such that the general solution 
is 

u = R 1 {p) + p'R 2 (p), (169) 

in which p = p(£ — Co, 32, 33) is the Weierstrass elliptic function characterized by (|24|l . 

If the genus is zero, there exists a (possibly zero) constant a and a rational function R such 
that the general solution is 

u = i?(e Q «), (170) 
with the degeneracy u = i?(£) in case a is zero. 
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